Observation of intermittency in wave turbulence 
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We report the observation of intermittency in gravity-capillary wave turbulence on the surface of 
mercury. We measure the temporal fluctuations of surface wave amplitude at a given location. We 
show that the shape of the probability density function of the local slope increments of the surface 
waves strongly changes across the time scales. The related structure functions and the flatness are 
found to be power laws of the time scale on more than one decade. The exponents of these power 
laws increase nonlinearly with the order of the structure function. All these observations show the 
intermittent nature of the increments of the local slope in wave turbulence. We discuss the possible 
origin of this intermittency. 

PACS numbers: 47.35.-i, 47.52. +j, 05.45.-a 



One of the most striking feature of turbulence is the 
occurrence of bursts of intense motion within more quies- 
cent fluid flow. This generates an intermittent behavior 
P, One of the quantitative characterization of in- 
termittency is given by the probability density function 
(PDF) of the velocity increments between two points sep- 
arated by a distance r. Starting from a roughly Gaus- 
sian PDF at integral scale, the PDFs undergo a contin- 
uous deformation when r is decreased within the inertial 
range and develop more and more stretched exponen- 
tial tails Deviation from the Gaussian shape can be 
quantified by the flatness of the PDF. The origin of non- 
Gaussian statistics in three dimensional hydrodynamic 
turbulence has been ascribed to the formation of strong 
vortices since the early work of Batchelor and Townsend 
[l|. However, the physical mechanism of intermittency is 
still an open question that motivates a lot of studies in 
three dimensional turbulence ^J] . Intermittency has been 
also observed in a lot of problems involving transport by 
a turbulent flow for which the analytical description of 
the anomalous scaling laws can be obtained [5|- 

It has been known since the work of Zakharov and col- 
laborators that weakly interacting nonlinear waves can 
also display Kolmogorov type spectra related to an en- 
ergy flux cascading from large to small scales ji', . These 
spectra have been analytically computed using perturba- 
tion techniques, but can also be obtained by dimensional 
analysis using Kolmogorov-type arguments f8l|. More re- 
cently, it has been proposed that intermittency correc- 
tions should be also taken into account in wave turbu- 
lence [§| and may be connected to singularities or co- 
herent structures [1, [13] such as wave breaking jlTi] or 
whitecaps [1] in the case of surface waves. However, in- 
termittency in wave turbulence is often related to non 
Gaussian statistics of low wave number Fourier ampli- 
tudes [lO|] , thus it is not obviously related to small scale 
intermittency of hydrodynamic turbulence. Surprisingly, 



there exist only a small number of experimental stud- 
ies on wave turbulence 12, l^, 13, 15, 3] compared to 



hydrodynamic turbulence, and to the best of our knowl- 
edge, no experimental observation of intermittency has 
been reported in wave turbulence. 

In this letter, we report the observation of an intermit- 
tent behavior for gravity-capillary waves on the surface 
of a layer of mercury. We show that we need to compute 
the second-order differences of the surface wave ampli- 
tude in order to display intermittency. We observe that 
the shape of their probability density function changes 
strongly across the time scales (from a Gaussian at large 
scales to a stretched exponential shape at short scales). 
This short-scale intermittency is conflrmed by comput- 
ing the structure functions for various time scales. The 
structure functions of order p (from 1 to 6) and the flat- 
ness are found to be power laws of the time scale on more 
than one decade. The exponents of the power laws of the 
structure functions are found to depend nonlinearly on 
p. All these observations show the intermittent nature of 
the local slope increments of the turbulent surface waves. 

The experimental setup has been already described 
elsewhere It consists of a square vessel, 20 x 20 cm^, 
filled with mercury up to a height of 18 mm. Mercury is 
chosen because of its low kinematic viscosity (one order 
of magnitude smaller than that of water), thus reducing 
wave dissipation. Note however that similar qualitative 
results to the ones reported here are found when chang- 
ing mercury by water. Surface waves are generated by 
the horizontal motion of one rectangular (13 x 3.5 cm^) 
plunging Plexiglas wave maker driven by an electromag- 
netic vibration exciter. The wave maker is driven with 
random noise excitation, supplied by a function genera- 
tor, and selected in a frequency range 0-6 Hz by a low- 
pass filter. The rms value of the velocity fluctuations 
of the wave maker is proportional to the driving voltage 
Urms applied to the vibration exciter. Surface waves are 
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generated 3 cm inward from one vessel wall. The local 
vertical displacement of the fluid is measured, 7 cm away 
from the wave maker, by a capacitive sensor. The sensor 
allows wave height measurements from 10 /xm up to 2 cm 
with a 20 mm/V sensitivity and a 0.1 ms response time. 

A typical recording of the surface wave amplitude, "^{1), 
at a given location is displayed in the inset of Fig. [T] as 
a function of time. The surface strongly fluctuates with 
a large distribution of amplitudes (see afterwards). The 
mean value of the amplitude is close to zero. In order to 
characterize the statistical properties of such a signal (in- 
set of Fig. [T]), r]{t) is recorded by means of an acquisition 
card with a 1 kHz sampling rate during 3000 s, leading 
to 3 X 10^ points recorded. The power spectrum and the 
probability density function (PDF) are then computed. 

At high enough forcing, the signature of a wave turbu- 
lence regime is observed 17[: a scale invariant spectrum 
with two power-law frequency dependences (see Fig. [T]) 
and an asymmetric PDF (see Fig. [2]). The low frequency 
spectrum part ~ f~^'^ corresponds to the gravity regime, 
whereas the high frequency one /~^'^ corresponds to 
the capillary regime. For the present characteristics of 
the forcing, both power-law exponents are in fair agree- 
ment with weak turbulence theory predicting a power 
spectrum of the wave amplitudes ^ for gravity waves 
y, and ^ /~^^^^ for capillary waves [3]. However, the 
scaling has been also ascribed to cusps In ad- 
dition, as emphasized in [13], only the capillary regime 
is robust, the exponent for the gravity regime being 
strongly dependent of the characteristics of the forcing. 
The cross-over near 30 Hz corresponds to the transition 
between gravity and capillary wave turbulence spectra. 
At still higher frequencies (greater than 150 Hz), viscous 
dissipation dominates and ends the cascade of energy in- 
jected from large scale forcing. 

The statistical distribution of wave height, 77, at a given 
location is displayed in Fig. [51 At high enough forcing, 
the PDF is no longer Gaussian, and becomes asymmetric. 
The positive rare events such as high crest waves are 
more probable than deep trough waves. This also can be 
directly observed on the temporal signal ri{t) in Fig. [TJ 

To test the intermittent properties of a stochastic sta- 
tionary signal r](t), one generally computes the incre- 
ments 5rj{T) = 7]{t + t) — rj{t). The structure functions 
of the signal Sp(r) = {\5ri(T)\v) ^ (h(i + r) - ry(t)|P) are 
also computed to seek a possible scaling behavior with 
the time lag r Q, (•) denoting a temporal average. How- 
ever, if a signal has a steep power spectrum E{f) ~ /~" 
with n > 3 (e.g. in Fig. [1]), the signal is then at least one 
time differentiable, and the increments are thus poorly 
informative (since they are dominated by the differen- 
tiable component), and S2(t) = {[rjit + r) 



V{t)?) - 



whatever n 0, 113]. To test intermittency properties 
of such a signal, a more pertinent statistical estimator 
is related to the second-order differences of the signal, 
Ary(r) = ■q{t + t) - 2'q{t) + r]{t - t) The struc- 
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FIG. 1: Power spectrum of surface wave height, ri{t). Dashed 
lines have slopes -4.3 and -2.8. Inset: Typical recording of rj{t) 
at a given location during 10 s. {77) ~ 0. Forcing amplitude 
Urms = 0.4 V. Forcing frequency band < / < 6 Hz. 




FIG. 2: Probability density function of the normalized wave 
height, ri{t)/arf. Standard deviation (t,, = ^^irf) — 2.6 mm, 
flatness {ri*)/{'iff = 4, and skewness {Ti"^) /{-rff^ = 0.65. 
Gaussian fit with zero mean and unit standard deviation 
( ). Same forcing parameters as in Fig. [1] 



ture functions are then defined as Sp{T) = {\/S.t](t)\p) . 
Note that other more complex estimators exist based on 
wavelet analysis or on inverse statistics |23|] of smooth 
signals. 

The probability density functions of the second-order 
differences of the surface wave height, A77(r), normalized 
to their respective standard deviation Ur, are plotted in 
Fig.[3]for different time lags 6 < r < 100 ms, the correla- 
tion time of ri{t) being Tc 63 ms. A shape deformation 
of the PDFs of A77(t) /cr is observed with the time lag 
T. The PDF is nearly Gaussian at large t. When r is de- 
creased from this integral scale, the PDF's shape changes 
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continuously, and strongly differs from a Gaussian (see 
the PDF's tails in Fig. [3]). This is a direct signature of 
intermittency. The extreme fluctuating events (large val- 
ues of A'r]{T)/ar) are all the more likely when the time 
scale T is short. Thus, the signal of the surf ace- wave 
amplitude displays intermittent bursts during which the 
slope varies in an abrupt way within a short time. The 
second-order differences of the wave-amplitude signal are 
indeed related to intermittency of the local slope incre- 
ments of the surface waves. 




FIG. 3: Probability density functions of second-order differ- 
ences of the wave height [ri{t+T)—2ri{t)+ri{t—T)]/aT for differ- 
ent time lags 6 < r < 100 ms (from top to bottom). Gaussian 
fit with zero mean and unit standard deviation ( ). Cor- 
relation time Tc — 63 ms. Each curve has been shifted for 
clarity. Same forcing parameters as in Fig. [1] 

Figure H] shows the structure functions Sp{T) of the 
second-order differences of the wave amplitude as a func- 
tion of the time lag r. For 5 < r < 50 ms, all the 
structure functions of order p = 1 to 6 are found to be 
power laws of r, 5p(r) r^p, where is an increasing 
function of the order p. When t > Tc, Spi^r) is found to 
saturate (e.g. towards 2(77^), for p=2) as usual ^] (data 
not shown). To quantify the intermittency of the sig- 
nal (i.e. the PDF shape deformation across the temporal 
scales), the dependence of the flatness, S4/S2, as a func- 
tion of T is displayed in the inset of the Fig. [H At large 
T, the flatness is close to 3 (the value for a Gaussian) and 
increases up to 26 at the shortest r, corresponding to a 
much flatter PDF (see Fig. [3]). The flatness is a power 
law of the time scale: S^/S^ ~ t" with c = -0.88 ± 0.03. 

The evolution of the exponents of the structure 
functions as a function of p is shown in Fig. [5] from the 
slopes of the log-log curves in Fig. |4l ^p is found to be 
a nonlinear function of p such that ^p — cip — ^p^ with 
ci = 1.65 ± 0.05 and C2 = 0.2 ± 0.02. The value of ci is 
related to the exponent of the low frequency spectrum, 
— 2ci — 1 « —4.3 ± 0.1. As said above, this value of ci 



can be related to the cusps observed on the fluid surface. 
They correspond to discontinuities in the vertical veloc- 
ity V of the surface, thus leading to /^^ spectrum, i.e., 
{[v{t + t) — v{t)]'^) cx T. This leads to the dimensional 
estimate (\r](t + T)-2ri{t) + r](t-T)\P) cx tP+p/^ = r^P/^, 
in fair agreement with the measurements for p = 1 and 
2 (see Fig. [5]). The nonhnearity of ^p (cn ^ 0) is an- 
other direct signature of intermittency [3|. This inter- 
mittency is observed for 20 < 1/t < 200 Hz, that is, for 
the capillary wave regime. The so-called intermittency 
coefficient C2 can be also deduced from the measurement 
of the flatness as a function of t (inset of Fig. [4|). In- 
deed, inserting the expression of into F = S4/S2 with 
Sp{t) ^ T^p, leads to F ~ r~^'^^. Thus, our measure- 
ments give C2 = — c/4 = 0.22 ± 0.008 which is in agree- 
ment with the value of C2 deduced from the exponents 
of the structure functions. Note that this intermittency 
coefficient is robust when using a third-order increment 
processing (see <>-symbols in Fig.[5|). 
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FIG. 4: Structure functions <Sp(r) of the second-order differ- 
ences of the wave amplitude as functions of the time lag r, 
for 1 < p < 6 (as labeled). (— ): Power law fits, iSp(r) ~ r^^, 
where the slopes depend on the order p (see Fig. [5]). Curves 
has been shifted for clarity. Inset: Flatness S4/S2 as a func- 
tion of r. ( — ): Power law fit with a slope —0.88. Correlation 
time Tc — 63 ms. Same forcing parameters as in Fig. [T] 

Figure \5\ shows also the exponents of the structure 
functions computed from the first-order differences of the 
signal, (|77(t + t) — ri{t)\P) ^ r^f, (open circles in Fig. 
[5]). One can thus compare with theoretical predictions 
of weak turbulence 0| (solid circles in Fig. [5]). Dimen- 
sional analysis for weak capillary wave turbulence gives 
^p — llp/12 (solid line). Although the experimentally 
measured slope is slightly smaller than 11/12 ~ 0.92, it 
is too close to 1 in order to display intermittency just by 
computing the structure functions from the signal incre- 
ments. Indeed, as said above, since our signal spectrum 
is very steep, the signal increments are poorly informa- 
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FIG. 5: Exponents of the structure functions as a function 
of p. 5p computed from the (□) second-order differences (from 
the slopes of Fig. or from the (0) third-order differences, 

and fitted by ( ) £,p — cip — with ci = 1.65 and C2 = 

0.2. ( — ) Theoretical prediction using dimensional analysis 
5p = 3p/2. (o) computed from the first-order increments, 
and fitted by (— •) (,p = 0.85p. (•) Theoretical points (from 
the first-order increments) [24| and dimensional estimate ^p — 
llp/12 (-). 



tive, and the exponents of the structure functions of 
the signal increments, {\r]{t + r ) — t]U)\p) ^ t^" , are ex- 
pected to be such that ip ^ p |19l . |20l| . This shows that 
second-order differences should be computed in order to 
test intermittency properties. 

We have reported that short-scale intermittency oc- 
curs on the second-order differences of surface wave am- 
phtude. As previously proposed, intermittency could be 
related to coherent structures on the fluid surface [lo| . 
such as wave breaking or whitecaps 8]. Here, wave 



breaking or whitecaps do not occur, but cusps are ob- 
served on the fluid surface. However, we do not presently 
have a theory that determines for large p. We think 
that the observation of small scale intermittency in our 
system that strongly differs from high Reynolds number 
hydrodynamic turbulence is of primary interest. It can 
indeed motivate explanations of intermittency different 
than the ones considering the dynamics of the Navier- 
Stokes equation or the existence of coherent structures. 
A more general explanation can be related to the proper- 
ties of the fluctuations of the energy flux that are shared 
by different systems displaying an energy cascade. 
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